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Geometric frustration in quantum magnetism refers to that magnetic interactions on different 
bonds cannot be simultaneously minimized. Usual Cooper pairing systems favor uniform spatial 
distributions of pairing phases among different lattice sites without frustration. In contrast, we 
propose "frustrated Cooper pairing" in non-bipartite lattices which leads to the supersolid states of 
Cooper pairs. Not only the amplitudes of the pairing order parameter but also its signs vary from 
site to site. This exotic pairing state naturally occurs in the p-orbital bands in optical lattices with 
ultra-cold spinless fermions. In the triangular lattice, it exhibits an unconventional supersolid state 
with the /-wave symmetry. 
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I. INTRODUCTION 

Frustration is one of the fundamental challenges in 
classic and quantum magnetism 1 . For the antiferromag- 
netic states in non-bipartite lattices, such as triangular, 
Kagome and pyrochlore, it is impossible to simultane- 
ously minimize the magnetic energy of each bond. Con- 
sequentially, the ground state configurations are heavily 
degenerate. The enhanced thermal and quantum fluctu- 
ations strongly suppress spin ordering. The ultimate or- 
derings often occur at much lower temperatures than the 
energy scale of the antiferromagnetic coupling through 
the "order from disorder mechanism" 2,3 . Furthermore, 
frustration provides a promising way to reach the spin 
liquid states, which exhibit exotic properties including 
topological ordering and fractionalizationi£. 

In the usual superfluid states of paired fermions and 
bosons, a uniform distribution of the superfluid phase 
over the lattice sites is favored in order to maximally 
facilitate phase coherence. However, frustration can 
indeed occur under certain conditions. It has been 
found that in the disordered superconductors near the 
superconductor-insulator transition, the fluctuations of 
the superfluid density can result in frustrated Joseph- 
son coupling among superconducting grains 6 - 8 . Recently, 
the striped superconductivity^^ has been proposed for 
the high T c compound La2_ a; Ba 2 ;Cu04. The Josephson 
coupling between two adjacent superconducting stripes 
is like in the 7r-junction leading to the opposite signs 
of the pairing phases across the junction. The mecha- 
nism for the frustrated coupling arises from the inter- 
play between superconductivity and antiferromagnetism 
in doped Mott insulators. However, an intuitive pic- 
ture of the microscopic origin of this exotic phase is still 
needed. 

On the other hand, cold atom optical lattices have 
opened up a new opportunity to investigate novel fea- 
tures of orbital physics which do not exhibit in usual 
orbital systems of transition metal oxides. Bosons have 
been pumped into the excited p-orbital bands experimen- 
tally with a long life timeii~— . This metastable excited 



state of bosons does not obey the "no-node" theory and 
exhibits the unconventional superfluidity with complex- 
valued many-body wavefunctions breaking time-reversal 
symmetry spontaneously^—, which has already been 
observe d 12 ' 13 . For orbital fermions, large progress has 
been made in the p^-orbital bands in the hexagonal 
lattices, whose physics is fundamentally different from 
that in the p z -orbital system of graphene. The inter- 
esting physics includes the flat band structure^, the 
consequential non-perturbative strong correlation effects 
(e.g. Wigner crystal** and ferromagnetism 20 ), frus- 
trated orbital exchange interactional, quantum anoma- 
lous Hall effect 22 , and the unconventional /-wave Cooper 
pairings*. 

We are interested in bridging the above important 
research directions together by introducing frustration 
to Cooper pairing as a new feature of orbital physics. 
In this article, we propose the "frustrated Cooper pair- 
ing" in the p^^-band of the non-bipartite optical lattices 
with spinless fermions. Due to the odd parity nature 
of the p^y-orbitals, the Josephson coupling of the on- 
site Cooper pairing is frustrated. In the strong coupling 
limit, the super-exchange interaction of the pseudo-spin 
algebra composed of the pairing and density operators is 
described by the "antiferromagnetic" Heisenberg model 
with the Ising anisotropy. It results in the coexistence 
of charge density wave and superfluidity of Cooper pairs 
with a non-uniform phase pattern. This supersolid state 
of Cooper pairs exhibits the /-wave pairing symmetry 
in the triangular lattice within a large range of particle 
density. 

Before we move on, let us explain some conceptual sub- 
tleties. One might wonder how to justify the validity of 
"frustration" of Cooper pairs which are usually an ex- 
tended objects. Indeed, frustration is mostly commonly 
defined in antiferromagnetism of local spin moments. 
However, frustration does not necessarily mean "on-site" 
physics even in the context of antiferromagnetism in non- 
bipartite lattices. For example, antiferromagnetic orders 
can be considered as pairing between particles and holes 
in the spin triplet channel carrying nonzero momentum, 
i.e., spin density waves. In the strong coupling limit, the 
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particle-hole pairs are strongly bound to be on-site, then 
the physics reduces to local moments described by the 
Heisenberg model. However, in the weak and intermedi- 
ate coupling regimes with small charge gaps, the systems 
are still locally itinerant. The spatial extensions of the 
particle-hole bound states are beyond one lattice site. If 
the lattice is non-bipartite, we still have frustrated mag- 
netism with extended particle-hole pairs. For example, 
this picture applies to the intermediate coupling regime 
of the Hubbard model at half-filling in the triangular lat- 
tice. In our case, we will consider the Cooper pairing at 
intermediate and strong coupling regimes. In the strong 
coupling limit, Cooper pairs are bound on a single site, 
whose exchange physics can be described by the antifcr- 
romagnetic pseudospin Heisenberg model in the charge 
channel. In the intermediate coupling regime, although 
a Cooper pair is an extended objects covering several 
lattice constants, its location can still be defined by its 
center of mass. The associated physical quantity is the 
pairing order parameter at each lattice site in the mean- 
filed theory. This physics can be best explained in terms 
of the anomalous Green's function F(R,r;uj), where R 
is the center of mass coordinate, and r is the relative 
coordinate. The order parameter A(i?) corresponds to 
F(R, 0; to = 0), while the size of Cooper pairing is deter- 
mined by the decay length of F(R, f; ui = 0) with respect 
to r. In our context, frustration refers to center of mass 
motion of A(R). 

This paper is organized as follows. The model Hamil- 
tonian and the band structure are introduced in Sect. [II] 
The strong coupling analysis is given in Sect. IIII1 The 
mean-field theory analysis is presented in Sect. IIV1 and 
the /-wave supersolid state is presented in [V] Conclu- 
sions are given in Sect. [VT] 



II. THE MODEL HAMILTONIAN 




FIG. 1: (a) The <r-bonding and 7r-bonding of the p-orbitals 
have opposite signs due its odd parity nature. This gives 
rise to the ' 'antiferromagnetic"-like exchange in the change 
channel with attractive interactions as expressed in Eq. [5] 
which is frustrated in the triangular lattice, (b) The first 
Brillouin zone of the triangular lattice is a regular hexagon, 
if i,2 = (±|j,0) represent two non- equivalent vertices, and 
other four are equivalent to Ki,2- 

We take the 2D triangular lattice as an example, which 
has been constructed experimentally by three coplanar 



laser beams^. The optical potential on each site is ap- 
proximated by a 3D anisotropic harmonic potential with 
frequencies uj z 3> lo x = ui v . After the lowest s-band is 
fulfilled, the active orbital bands become p x ,y The p z - 
band remains empty and is neglected. The free part of 
the pa-^-orbital band Hamiltonian in the triangular lat- 
tice filled with spinless fermions reads 

Ha = t \\ {PLfjPLf+aeij + h.cA 

r ,1— 1~3,£T 

- t± ^ (pT,f,iPT,P+ae z ,i + h.C.J 

- fx^n^, (1) 

where r runs over all the sites; €\ — e x , £2,3 = — ie x ± 

■^e y are the three unit vectors along bond directions. 
PL.i = (p x e x + p y e y ) • Bi are the longitudinal projections 
of the p-orbitals along the ii direction. More explic- 
itly, p L<1 = p x and Pl, 2(3) = ~\Px ± ^rPy The trans- 
verse projections of the p-orbitals along the bond read 
PT,i = {p x e x +P v By) ■ (z x St), n? = p\p x + p\p y is the 
particle number operator; \x is the chemical potential. 
The cr-bonding t\\ and 7r-bonding t± describe the hop- 
pings between p-orbitals along and perpendicular to the 
bond direction, respectively, as depicted in Fig. [T]A. tt\ 
is positive due to the odd parity nature of the p-orbitals, 
which is scaled 1 below. t± is usually much smaller than 
tii because of the anisotropy of the p-orbitals. The first 
Brillouin zone (BZ) of the triangular lattice is a regular 
hexagon as plotted in Fig. [T]B. The edge length of the 
first BZ is 4^, where a is the lattice constant. 

3a 5 




FIG. 2: The Band structure of the Hamiltonian Eq. [T] Two 
bands touch at Ki t z with the Dirac spectra, and at the center 
of the BZ with the quadratic spectra. 

The band structure of the noninteracting Hamiltonian 
Eq. [TJ is displayed in Fig. [5] with the value of t±/tu 
chosen as 0.2. The band width is around 6ty. There 
is no particle-hole symmetry with respect to the zero- 
energy point, which hints an asymmetric phase diagram 
with respect to half-filling for the interacting Hamiltonian 
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introduced in Sect. Mil In momentum space, we define 
the two-component spinor ijj(k) = (p x (k),p y (k)) T for the 
p x and Py-orbitals. The Hamiltonian Eq. Q] becomes 

H = ^&(%){Hc0$) - iri a p}4f}(k), (2) 
k 

where the matrix kernel H a p(k) takes the structure of 
H(k) = /(£) + <7i(£)n+53(£)T 3 , (3) 

where n,3 are the Pauli matrices defined for the basis 
of p x , p y for the spinor of ip(k); the expressions of f(k), 
gi(k) and gi (k) arc 

f(k) = (1 - <_L)^cosfc • ei, 

i 

gi(k) = - — (1 + tj_)(cosfc ■ e 2 - cosfc • e 3 ), 

g 2 (k) = -^(1 + tj_)(cosfc ■ e 2 + cosfc ■ e 3 ), 

+ (1 +t±)cosk-e 1 . (4) 

The diagonalization of H (fc) gives rise the dispersions of 
two bands as 

E ± = f(k)±(l + t ± ) 

x /^cos 2 fc-ei— cosfc • e a cosfc • e&. (5) 

V i l<a<b<3 

These two bands touch each other at with the Dirac 
cone-like spectra, and at the center of the BZ with the 
quadratic spectra. When the Fermi energy is located 
at the Dirac points, the other band contributes a large 
connected branch Fermi surface, thus its contributions to 
thermodynamic quantities dominate over those from the 
Dirac points. 

The topology of the Fermi surfaces varies at different 
filling levels. The energy minima of Eq. [T]are three-fold 
degenerate located at the middle points of the BZ edges. 
The middle points of the opposite edges are equivalent 
up to a reciprocal lattice vector. Around the band bot- 
tom, the Fermi surfaces only cut the first band and form 
three disconnected elliptical pockets. As filling increases, 
these pockets become connected forming a large Fermi 
surface around the center of the BZ. At the same time, 
the two Dirac cones contribute two Fermi surfaces around 
the K\ 2 points, which shrink to two points when Fermi 
energy is right at the Dirac points. As approaching the 
band top where a quadratic band touching exists, there 
are two Fermi surfaces around the center of the BZ. 

III. STRONG COUPLING ANALYSIS 

It is well-known that in the Mott-insulating state of 
the positive- U Hubbard model at half- filling, its low en- 
ergy physics lie in the magnetic channel, which is cap- 
tured by the antiferromagnetic Heisenberg models. In 




FIG. 3: The pattern of {G\fJ\G} at h = 0: the unit cell of 
three sites exhibiting the supersolid ordering. 

non-bipartite lattices (e.g. the triangular lattice), the an- 
tiferromagnetic exchange cannot be simultaneously min- 
imized for every bond, which leads to frustration. Sim- 
ilarly, for the negative-C/ Hubbard model, in the strong 
coupling limit, the low energy physics is described by the 
exchange interaction in the charge channel. On each site, 
the low energy states are the doubly occupied state and 
the empty state, which can be considered the pseudospin 
"up" and "down" state, respectively^. The Josephson 
coupling between neighboring sites plays the role of the 
ferromagnetic coupling in the xy-direction of the pseu- 
dospin, which favors a uniform phase distribution. In 
other words, the exchange interaction in the charge chan- 
nel is unfrustrated. 

This situation is fundamentally changed for the Hub- 
bard model of spinless fermions in p^-orbitals based on 
the band structure in Eq. [TJ We add the attractive Hub- 
bard interaction between spinless fermions in the p x . y - 
orbital bands 

Hint = -U / J'V,p«>V,i>„) ( 6 ) 

r 

where U is positive. The frustrated nature of Cooper 
pairing can be easily explained in the strong coupling 
limit of U ^> in. Similarly to the usual negative U 
Hubbard model, we construct the pseudospin algebra 
denoted^ as 

Vx = ^(ptPy+PyPx), Vy = -^(plPy-PyPx), 

Vz = ^(np-l)- (7) 

Up to a normalization factor, they are the real and imag- 
inary parts of the pairing operator, and the particle den- 
sity operator, respectively. The low energy Hilbert space 
in each site consists of the doubly occupied state and the 
empty state, which are eigenstate of rj z with eigenvalues 
±i, respectively. 

This super-exchange interaction of the pseudo-spin has 
a remarkable feature that a 7r-phase difference is favored 
between pairing order parameters rj x _ y on neighboring 
sites. As a pair hops to the neighboring site, it gains a 7r- 
phase shift, because the a-bonding and 7r-bonding terms 
are with opposite signs as depicted in Fig. [T]A. The per- 
turbation theory gives rise to the anisotropic "antiferro- 
magnetic" Heisenberg model (AAFHM) in the "external 
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FIG. 4: The average fermion number per site n versus the 
chemical potential /i at U/tii = 6 and t±/t\\ = 0.2. 



magnetic field" , 

He ff = E J ^,v{Vx{i)Vx(j) + Vy(i)Vv(J)} + JzVz{i)Vz(j) 

ij 

- ft£>(0, (8) 

i 

where h = 2/j is the external magnetic field; the exchange 
read 



constants J x ^ y and J. 



Jx,y — 



44 i in 



U 



Jz = 



2(4l + 4 2 ) 

c7 ' 



(9) 



The Ising anisotropy in Eq. [3] is because J z > J x ,y Sim- 
ilar models apply to pairing problem of spinless fermions 
in the p x ijr orbital of the bipartite lattices of square 2 ^ and 
hexagonal 2 ^ which are not frustrated because a canonical 
transformation can change J x . y to — J x _ y . 

Eq. [5] can be interpreted as a hard core boson model 
with the frustrated hopping J XtV and the nearest neighbor 
repulsion J z . It has been studied at the zero external field 
in Ref j 29 ' 30 which shows a supersolid orderin g 31 ' 32 with 
a three-site unit cell as depicted in Fig. [3j The competi- 
tion between charge density wave and supersolid ordering 
in optical lattices has also been studied in Refs. |33ll3 4| . 
Site A has no superfluid component, i.e., 77 || z; sites B 
and C develop superfluid orders with a 7r-phase differ- 
ence. However, the experimental realization of hard core 
bosons with frustrated hopping is difficult. In compari- 
son, our idea of the frustrated Cooper pairing of fermions 
is very natural in the p X)l/ -orbital bands. Furthermore, 
previous studie a 29 ' 30 focus on the pseudo-spin model Eq. 
Incompletely neglecting the fermion degree of freedom. In 
the following, instead of using Eq. [SJ we directly study 
the Cooper pairing problem with the fermion Hamilto- 
nian in the entire filling range from to 2. 



IV. MEAN-FIELD THEORY AT 
INTERMEDIATE COUPLINGS 

Below we will focus on the intermediate coupling 
regime and perform the self-consistent mean-field the- 
ory to the fermionic Hubbard model of Eq. [T] and Eq. 



[6] Unlike the positive-?/ Hubbard model with doping, in 
which the mean-field theory is unreliable, in our case of 
negative-?/ Hubbard model, the mean-field theory gives 
qualitatively correct results for the competition between 
Cooper pairing and charge density wave (CDW). For a 
detailed review, pleas refer to Ref. HE 

To decouple Hi nt , we assume the pairing and CDW 
ordering taking an enlarge unit cell of three sites, and 
define 

A/ = {G\pf> e i, v p? e i, x \G), Nr = ^{G\h reI \G), (10) 

where I = A, B,C refer to the sublattice index; (G\...\G) 
means the average over the mean-field ground state. The 
mean-field interaction Hamiltonian becomes: 



H mf = -U 



E 

?J=A,B,C 



{A*lPr£LyPr£l,x + h.c} 



N l{Pr£l,xPreI,x + pl ehy PreI.,y} 



(11) 



Combining Eqs. [Tl and [TT1 and performing Fourier trans- 
formation to momentum space, we can obtain the result- 
ing mean-field Hamiltonian as: 



H s (k) D(k) 
t(fc) -H* s (-k) 



*(*0 (12) 



where y]p means the summation only cover half of the 
reduced Brillouin zone; $(k) is defined as 



(13) 



where <f>(k) = [p A ,x{k),PA, v (k) j PB.x 

(k) (k) 

•> PC,x 

Pc,y(k)]', H s contains the free Hamiltonian Eq. Q] com- 
bined with the CDW decoupling; D(k) is the pair- 
ing part. The order parameters are obtained self- 
consistently. The above definition of order parameters 
are related to the pseudospin operators through 

(G\ Vx (fe I)\G) = Re A/, (G\r) y (r e I)\G) = ImA,, 
(G\r) z (r£l)\G) - Nj — 1/2. (14) 

Different from the ordinary BCS problem, the pairing 
of Eq. [T] is not an infinitesimal instability but occurs at 
the finite attraction strength. It is because the eigen- 
states of the two time-reversal partners with momentum 
k and —A: of the free Hamiltonian Eq. [T]have the same 
real polar orbital configuration. This suppresses pairing 
at weak interactions because attraction only exits in or- 
thogonal orbitals. With intermediate and strong interac- 
tions, pairing can occur between different bands. Below 
we present results for tx/tu = 0.2 and an intermediate 
coupling and U/tn = 6. This corresponds to the effective 
AAFHM with the Ising anisotropy of J z /J x ,y = 2.6. 

We discuss our mean-field results in terms of the pseu- 
dospin orientations at the three sublattices A, B, C. Fig. 
Ushows the total fermion number per site n = (nA + ns + 
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FIG. 5: The real space configurations of pseudospin ff on 
the xz plane at various fillings n from (A) to (I). A) and (I) 
indicate fully polarized states. (B) and (H) show three titled 
vectors, where two of them have a relative 7r-phase to the third 
one. (C) and (G) depict the CDW insulating state. (D) and 
(F) exhibit an umbrella-like shape with opposite orientation. 
Two of them have a 7r-phase difference and the third one does 
not own superfluid component. (E) denotes an intermediate 
configuration between (D) and (F). 



(A) n=0.68 (B) n=0.7 

I L ! 

(D) n=0.75 (E) n=0.775 



(C) n=0.725 




(M) n=1 



FIG. 6: The smooth evolution of the two opposite "umbrella"- 
like configurations from n = 0.68 to n = 1. 



nc)/3 as a function of chemical potential /i, which is the 
counterpart of the magnetization in the AAFHM. The 
first prominent feature is the plateaus occurring at n = | 
and |, which is corresponding to those at (G\r) z \G) = ± ^ 
observed in the study of classical ground state of the 
AAFHM. These two plateaus are corresponding to CDW 
insulating states without superfluidity. As shown in Fig. 
[5j the corresponding pseudospin orientation for CDW in- 
sulating states is that all the pseudospins are fully polar- 
ized along the z axis with two of sublattice along the 
same direction and the remaining one along the opposite 
direction. 

Although Fig. [4] resembles the behaviors of the mag- 
netization obtained by the AAFHM 3 —, two major differ- 
ences exist. First, the widths of the two CDW plateaus 
in Fig. 2] are different while those of the AAFHAM are 
the same. We attribute this discrepancy to the differ- 
ent symmetry properties between the AAFHM and the 
Hubbard model of Eq. [T] and Eq. [5] with the asymmet- 
ric band structure shown in Fig. [2J The AAFHM has 
the symmetry of the rotation of 180° around the x-axis, 
i.e., r/ x — > i] x , rj y ,r/ z —} —i] y ,—r] z and h — > —h. Such an 
operator corresponds to the particle-hole transformation 
at the fermion level as p x — > ip^ and p^ ip x , which 



is not kept in the triangular lattice. As a result, for the 
AAFHM, the magnetization should be an odd function 
with respect to h so that the lengths of the plateaus are 
the same. This kind of behavior is not expected in Fig. 
|U The other difference is that at h = 0, the ferrimagnetic 
state is found in the AAFHM, and our results show the 
"paramagnetic" behavior, i.e., there is no jump around 
n = 1. This is due to the quantum fluctuations arising 
from the singly occupied states as discussed below. 

Fig. [5] plots the pseudospin orientations on three sub- 
lattices at a series of filling levels. Except the CDW 
insulating states at n — | , | , we find that the pseu- 
dospins have non-zero (G\r] x \G} and (G\r] z \G) in the most 
part of the phase diagram, indicating the frustrated su- 
persolid states with non-uniform Cooper pairing density 
and phase. Moreover, the phase diagram can be well- 
understood by the rotations of pseudospin orientation un- 
der the magnetic field h — 2//. At n = 0, h is large along 
the —z direction so that all the pseudospins are com- 
pletely polarized. As n increases, the magnitude of h de- 
creases so that the pseudospins gradually rotate upward 
with one of the pseudospins (ff(A)) having much faster ro- 
tating rate. They become polarized along the z-direction 
as arriving at the CDW insulating state n = | with one 
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FIG. 7: The /-wave pairing pattern in (A) real space and 
(B) momentum space. The rotation of 60° around the center 
site O in A) (denoted by the hollow circle) in real space or 
around the center of the reduced BZ in momentum space is 
equivalent to reverse the sign of the pairing order parameters. 

pointing up and the other two pointing down. The mag- 
nitudes of T)a,b,c are smaller than \ due to quantum 
fluctuations. As n increases further, since rj z {A) can not 
increase anymore, rf(B) and 77(C) gradually turn upward 
leaving if a unchanged forming an umbrella configuration. 
Such a state is a coexistence of superfluidity and CDW, 
thus is a supersolid state. 

After a critical value n c ~ 0.7, all the pseudospins 
start to rotate simultaneously and continuously evolve 
between the two umbrella configurations with opposite 
orientations depicted in (D) and (F) in Fig. [5l respec- 
tively. A detailed process of evolution is plotted in Fig. 
[6] from n = 0.68 to n = 1. This continuous evolution 
of the ground state is not present in the AAFHM since 
its ground state is ferrimagnetic with non-zero magneti- 
zation at zero field, which corresponds to n ^ 1 in our 
model. This deviation is because the AAFHAM model 
is only justified at the strong coupling limit. The larger 
kinetic energy in this region leads to the less stringent 
assumption of the strong-coupling. Consequently, the 
quantum fluctuations arising from the singly occupied 
states are enhanced, which are in disfavor of CDW but 
in favor of uniform superfluidity The continuous evolu- 
tion also explains why there is no jump at n = 1 in Fig. U 
Finally, the rest part of the phase diagram can be easily 
understood by rotating all the pseudospins upward, and 
eventually all the pseudospins are fully polarized along 
+z direction at n = 2. 



V. THE /-WAVE SUPERSOLID STATE 

One remarkable feature of the frustrated Cooper pair- 
ing that we are studying is that it can give rises to an un- 
conventional type supersolid state exhibiting non-s-wave 
symmetry. For example, in Fig. [5] for a wide region of n 
(0.67 < n < 0.7, 1 < n < 1.3), we find (G\r) x (A)\G) = 
and (G\r) x (B)\G) = -{G\<q x {C)\G) = A. As shown in 
Fig. [7] A, the signs of the pairing order parameter are 
opposite in sublattices A and B. As a result, a spatial 




FIG. 8: A) The reduced Brillouin zone (RBZ) associated with 
the enlarged three-site unit cell compared with the original 
BZ. The six vertices of the RBZ are located at the centers 
of the six regular triangles composed of the center of the BZ 
and the vertices of the original BZ. B) The six bands in the 
reduced BZ with the three-site CDW pattern of ua = nc = 
1.52, and n B = 0.154 for t± = 0.2 and U = 6. The chemical 
potential fi is reset to zero. 

rotation around a site of sublattice C at 60° corresponds 
to flipping the sign of the order parameters, which indi- 
cates the /-wave pairing symmetry. 

The /-wave pairing symmetry is also manifest in the 
gap function structure in momentum space. We calculate 
the intra-band pairing functions A„„ in the momentum 
space by projecting the pairing potential in Eq. Qj]to the 
band eigen-basis as: 

/ 6 

£ E Km(k)MWm(-k)+h.C, (15) 
u m,n— 1 

where 

A nm (k) = [uHk)D(k)U*(-k)] nm , (16) 

U(k) is the unitary matrix such that 

U\k)H s (k)U(k)=diag[E 1 (k),...,E 6 (k)], (17) 

and H s (k), D(k) are given in Eq. [T2j We have confirmed 
that all six intra-band pairing functions have three nodal 
lines and sign changes under 60° rotation. 

As a specific example, we present the results of the 
self-consistent mean-field theory for the filling n = 1.07 
and other parameters t± = 0.2 and U — 6 as before. The 
system exhibits the three-site pattern of the CDW order 
as ua = nc = 1.52, and ub — 0.154, and the pairing or- 
der parameters in the real space are = —Ac = 0.345 
and Ab — 0. We plot the band structure with the above 
CDW order parameter but set the gap functions zero. 
The reduced BZ is only | of the original BZ, and there 
are six bands in total as plotted in Fig. [8] The chem- 
ical potential 71 is reset to 0, which lies in the gap be- 
tween the 4th and 5th bands and has no crossing with 
the band spectra. As a result, although the gap functions 
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have node lines due to the /-wave symmetry, the Bogoli- 
ubov excitations remains fully gapped. We plot the gap 
function of A44 in Fig. [7] B for demonstration purpose. 
The nodal lines are the three lines connecting the middle 
points of the opposite edges of BZ. Thus this is an uncon- 
ventional supersolid state of frustrated Cooper pairing 
with the /-wave pairing symmetry. Another interesting 
feature is that the gap function A44 even changes sign 
along the radial direction. 

It would also be instructive to compare our /-wave 
pairing supersolid state of Cooper pairs with the Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) state^Ml. Both 
cases exhibit non-uniform distributions of pairing phase 
in real space. However, the FFLO state completely 
breaks rotational symmetry. Its pairing pattern does not 
form a well-defined representation of the lattice point 
group in momentum space. In our case, it has a well- 
defined /-wave symmetry. 

We also consider the extreme anisotropy limit of the 
vanishing 7r-bonding strength, i.e., t± = 0. The bond 
superexchange only results in the J 2 -term at the second 
order perturbation level in Eq. [8] The leading order 
of the hopping of the Cooper pairs occurs through the 
three-site ring exchange 

Aff = -^J'[rix{i)Vx{j) + r 1 y{i)r lv {j)]r lz {k) (18) 

ijk 

where J' = § . The hopping is frustrated for a plaque- 
tte with only one site occupied, but it is unfrustrated for 
a plaquette with two sites occupied. This means that at 
low fillings the phase diagram does not change much from 



the case of nonzero t±, while the system finally evolves 
to a uniform pairing phase at n close 2. A more detailed 
analysis will be presented in a later publication. 



VI. CONCLUSION 

In summary, we introduce the concept of "frustrated 
Cooper pairing" of spinless fermions in the p-orbital band 
in optical lattices. The frustration occurs naturally from 
the odd parity of the p-orbitals and is a new feature of 
orbital physics. Exotic supersolid states of Cooper pairs 
with nonuniform distributions of pair density and phase 
are obtained with an unconventional /-wave symmetry. 
This opens up a new opportunity to study the physics of 
frustrated magnet by using the pseudo-spin algebra of the 
charge and pair degrees of freedom of Cooper pairs. This 
idea can also be applied to other even more frustrated 
lattices, such as Kagome and pyrochlore. In consider- 
ing the possibility of the existence of exciting spin liquid 
states therein, their counterparts in terms of "frustrated 
Cooper pairs" is another interesting direction for further 
exploration. 
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